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Abstract. The high frequency behaviour for random eigenfunctions 
of the spherical Laplacian has been recently the object of considerable 
interest, also because of strong motivations arising from Physics and 
Cosmology. In this paper, we are concerned with the high frequency 
behaviour of excursion sets; in particular, we establish a Uniform Central 
Limit Theorem for the empirical measure, i.e. the proportion of spherical 
surface where spherical Gaussian eigenfunctions lie below a level z. Our 
proofs borrows some techniques from the literature on stationary long 
memory processes; in particular, we expand the empirical measure into 
Hermite polynomials, and establish a uniform weak reduction principle, 
entailing that the asymptotic behaviour is asymptotically dominated by 
a single term in the expansion. As a result, we establish a functional 
central limit theorem; the limiting process is fully degenerate. 
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1. Introduction 

1.1. Background. Let A52 be the usual Laplacian on the sphere § 2 , i.e., 
given by 

1 9 r . „ a 1 id 2 

sin 06 y 06 J sir 9 oip z 

in the usual spherical coordinates (6, tp) , An orthonormal family of eigen- 
functions for the Laplacian is well-known to be given by the (complex- 
valued) spherical harmonics {¥jm} m =_j 1 > i- e - 

A s2 Y lm = -l(l+l)Y lm ,/ = 1,2,3,..., / Y lrn (p,(p)Y^(0,<p)wOMdv> = SZ,, 

Js 2 

where S^i is the usual Kronecker delta. For definiteness, we recall that in 
coordinates, the (complex- valued) spherical harmonics are defined by 



21 + 1 (I — m)\ 
Yim{0,y) = U 47r n + m y P im(cos8)exp(imip), 

jl+m 

P lm (x) = (-iHl-s^—^-l)', 
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{Plm(-)} denoting associated Legendre polynomials. The (21 + 1) spherical 
harmonics {Yj m } thus provide an orthonormal system for the space Hi of 
eigenfunctions corresponding to the eigenvalue —1(1 + 1), and we have the 
expansion 

oo 

L 2 (S 2 ) = 0^. 
i=i 

We shall equip each of the spaces {%i} with a rotation-invariant Gaussian 
probability measure, and focus on the asymptotic behaviour of the random 
eigenfunctions {/;(•)};=! 2 > i- e - 

l 

fl( x ) = ^2 a imyim(x) , X G S 2 , 
m=—l 

where the coefficients {a; m } are mean zero (complex) Gaussian random vari- 
ables such that Eai m = , Ea im ai m , = (21 + 1) _1 5™ , ai m = (-l) m a; _ m , 
whence 

Ef t (x) = , Eff(x) = (21 + l)E \a lm \ 2 = 1. 

Equivalently, // could be expanded in a real-valued basis; throughout this 
paper, however, we stick to the complex-valued representation, which has 
simpler symmetry properties. The random functions {/;(•)} are easily seen 
to be isotropic, i.e. for all g G SO(3) (the rotation group on M 3 ) we have 

fi(gx) = fi(x), 

= denoting identity in law, in the sense of stochastic processes. 

As well-known, the spherical Gaussian eigenfunctions emerge very nat- 
urally also in the analysis of random fields, since they provide the Fourier 
components in spectral representation expansions. The geometry of random 
fields has been of course extensively studied in the mathematical literature: 
we refer to the well-known book [1] for a complete list of references. Several 
statistics for Gaussian fields and their excursion sets have been investigated 
and many applications have also been implemented on experimental data, 
especially in a cosmological framework, see for instance jS] and the refer- 
ences therein. Geometric features of random spherical harmonics have been 
studied for instance by ([27], [25]) where sharp estimates for the variance of 
the nodal length (i.e. the boundary of excursion sets for z = 0) are provided. 
In this paper, we focus on a different characterization of the excursion sets, 
i.e. the behaviour of their empirical measure, to be defined below. 

1.2. Main results. Let us define the spherical harmonics empirical mea- 
sure as follows: for all z G (—00,00), 

(1) ^(z) := [ I(ft(x) < z)dx, 

Js 2 

where !(•) is, as usual, the indicator function which takes value one if the 
condition in the argument is satisfied, zero otherwise. In words, the function 
&l(z) provides the (random) measure of the set where the eigenfunction lie 
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below the value z. For example, the value of &i(z) at z = is related to the 
so-called Defect 

T>i := meas (/ ; _1 (0, oo)) — meas oo, 0)) 

by the straightforward transformation 

X>, = 47r-2$,(0). 

Of course, Air — &i{z) provides the area of the excursion set Ai(z) := 
{x : fi{x) > z} . The empirical measure can be viewed as a continuous ana- 
logue of the empirical distribution function for sequences of random vari- 
ables, and it is simply related to the first Minkowski functional from convex 
geometry, see p]. Clearly, for all z6i, 

E [*,(*)] =4tt*(s) j 

where <!>(.) is the cumulative distribution function of the standard Gauss- 
ian. Our aim in this paper is to study the distribution of {<&i(z) — &(z)} , 
uniformly w.r.t. z, asymptotically as I — > oo. The following lemma deals 
with the variance of &i(z) as / — > oo. 

Lemma 1.1. For every 

2 if A 2 1 , ^ f l °Sl 



Var($,(*)) = z l 4>{z) 2 - y + O z . /2 
where <j) is the standard Gaussian probability density function. 

In particular, for z / 0, Lemma 11.11 gives the asymptotic form of the 
variance as I — > oo. In contrast, for z = (this case corresponds to the 
Defect), this yields only a "o" -bound and one needs to work harder to obtain 
a precise estimate (see [29] ). Lemma 11.11 follows rather easily from the 
Hermite expansion approach in section [3l and Q, and we omit a formal 
proof. 

From Lemma 11.11 it is then natural to normalize &i(z) and define the 
spherical harmonics empirical process by 



(2) Gi(z) := Vl 



Hfl(x) < z)dx-{4ir x $(z)} 

s 2 



for I = 1,2,..., z G (—00,00). In a sense to be made rigorous later, we are 
then led to investigate the high-frequency ergodicity properties of Gaussian 
eigenfunctions (compare |17J). Indeed, as I diverges, ([2]) amounts to the 
distance between an average over a trajectory (f S 2^(fi(x) < z)dx) and an 
average over the ensemble determined by the probability law 

*(*) = / l(fi(x) < z)dP(u). 
Jn 

For G[(.), we establish a functional central limit theorem on a suitable Sko- 
rohod space. 

Instrumental for the main result is a uniform weak reduction principle, 
where we show that the behaviour of the empirical process is uniformly 
approximated by a quadratic polynomial in {//(•)} times a simple deter- 
ministic function. As a consequence, the weak limit is a fully degenerate 
random process, which can be realized as a single random variable times 
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the derivative of the Gaussian density. In the sequel, we use => to denote 
weak convergence in the Skorohod space D([— 00,00]), equipped with the 
sup-norm, and Goo (z) to label the mean zero, degenerate Gaussian process 

Goo(z) =z<j)(z)Z , Z~JV(0,1). 

Theorem 1.2. (The Uniform Central Limit Theorem) As I — > 00, 

Gi{z) => Goo(z). 

Considering the limiting expression for Gqo^), we remark the existence 
of a singularity for z = 0; indeed, here J g (0) = 0, whence it is easily seen 
that 

(3) G?,(0)-^ p 0. 

From ([3]), we immediately gather that the weak reduction argument no 
longer holds for the Defect (sec ([29 )): indeed, the limiting behaviour of 
the sequence {D{\ is considerably more complicated and will be considered 
in our forthcoming paper [29 1. 

The material in this paper attempts to provide a characterization as com- 
plete as possible on the high energy behaviour of the sample realizations 
of Gaussian eigenfunctions. A further natural question is to investigate 
whether the weak convergence result may enjoy some form of uniformity 
across different frequencies I. For the following result, => denotes weak con- 
vergence in the Skorohod space D {[—00, 00] x [0, 1]} , and Woo (2; r) denotes 
the mean zero Gaussian process with covariance function 

E[W O0 (zi;r 1 )W O0 (z 2 ;r 2 )] = {n A r 2 } z 1 z 2 (f>{z 1 )4>{z2) . 

Theorem 1.3. (Partial Sum Processes) As L — > 00 

W L (z;r) => Woo(z;r). 

Remark 1.4. It is possible to envisage some statistical applications of The- 
orem 11.31 for instance in the analysis of isotropic spherical random fields 
(such as those related to the analysis of the Cosmic Microwave Background 
radiation (CMB), see [7J for details). In such circumstances, the Gaussian 
eigenfunctions can be identified with the Fourier components of the field at 
frequency I, and it is straightforward to exploit Theorems \1.2\ and \1.3\ to con- 
struct tests for Gaussianity and isotropy, of Kolmogorov-Smirnov (Sl) or 
Cramer-von Mises (K£) type, by means of the statistics 

S L := sup sup I W L (z;r) | , K L := / / | W L (z; r)\ 2 dzdr 

z r Jo J-00 

where weak convergence ensures that 

Sl sup sup I Wqo (z; r ) I ,K L -t d I I |WooO;r)| dzdr, 

z r Jo J-00 



whence threshold values at given significance level can be easily tabulated. 
We refer for instance to |13| . |14j for more discussion on these issues. 
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1.3. Plan of the Paper. In Section 2 we discuss polynomial transforms of 
spherical random eigenfunctions and we establish their asymptotic behaviour 
under Gaussianity. Particular care is devoted to the derivation of exact 
asymptotic rates. In Section 3 we discuss the structure of our proofs and 
the significance of main results. In Section 4 and 5 we provide the proofs of 
Theorems (|1.2p and (j 1 . 3 [) . Some background material and technical lemmas 
are collected in an Appendix. 

1.4. Acknowledgements. We are grateful to Giovanni Peccati, Zeev Rud- 
nick and Mikhail Sodin for many stimulating discussions on these and related 
issues. This work was initiated during the workshop "Random Fields and 
Stochastic Geometry" held in February 2009, Banff, Canada, and the au- 
thors wish to use this opportunity to thank the organizers, Robert Adler 
and Jonathan Taylor, for the hospitality, and an excellent opportunity to 
exchange ideas and challenges. A substantial part of this research was done 
during the second author's visit to University of Rome "Tor Vergata" , and 
he would like to acknowledge the extremely friendly and stimulating envi- 
ronment in the institution, and the generous financial support. 

2. Polynomial transforms and their asymptotic behaviour 

Throughout this paper, we shall make extensive use of well-known re- 
sults on Hermite polynomials and higher-order moments of Gaussian ran- 
dom variables. To fix notation, we first recall the standard definition of 
Hermite polynomials 

the first few examples being provided by 

fli(aj) = x , H 2 (x) = x 2 - 1 , H 3 (x) = x 3 - 3x , H 4 (x) = x 4 - 6x 2 + 3 , 
H 5 (x) = x 5 - 10x 3 + 15a , H 6 (x) = x 6 - 15x 4 + 45x 2 + 15 , ... 
We recall also the basic formulae 

E[H q (Mx))]^0 

and 

E[H qi (Mx))H q2 (My))} = d% qi \ {EMx)My)r = 5% qi \I\((x,y))*, 
where 

Pi((*,v)) : =^7 E Y lm {x)Y lm {y)=E{f l {x)f l {y)} 

m=—l 

are the Legendre polynomials and (x,y) = cos(d(x,y)) where d(x,y) is the 
usual spherical distance on S 2 . For our arguments in the sequel, we shall need 
a few basic facts on the asymptotic behavior of the polynomial transforms 

(4) h l;q := [ H q {h{x))dx , for q = 1, 2, 3, .... 

Js 2 

Note first that Eh[. q = , for all I, q; also 

h\-\ = / fi(x)dx = , for all I > 1; 
Js 2 
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more generally, and hi- q = 0, when I and q are both odd, by symmetry 
conditions. The behaviour of quadratic transforms is discussed in [T7j, where 
non-Gaussian behaviour is also covered and uniform bounds on the speed of 
convergence to the limiting distribution are discussed, by means of a detailed 
analysis of the cumulants for spherical harmonics coefficients. We write here 
a different proof based on Legendre polynomials to introduce the techniques 
we shall use extensively, later in this paper. 

Proposition 2.1. (Proposition 2, |17j ) As I — > oo we have 

Vlhi i2 -+ d N(0,T). 



Proof. Note first that 



Var(h l]2 ) = Varlj 2 H 2 (fi( X ))dx\ = E U iW(s))cJx 
#[#2(/i0»0)#2(/i(y))] dxdy 



2 



S 2 xS 2 

Pf{t)dt 



1 2 







21 + 1 



In view of results [191 EQ] , to investigate the asymptotic behaviour of the total 
variation distance dTv{vlhi- t 2, Z) between Vlhi-2 and a standard Gaussian 
random variable Z ~ iV(0, 1), it is enough to show that the normalized 
fourth-order cumulants converge to zero, i.e., that 

cum A {hi. 2 ) 



{Varih^Y 

Using the well-known Diagram Formula (see e.g. [23J), we have 



cumi{hi- 2 ) = 48 / Pi((x,y))Pi({y,z))Pi({z,w))Pi({w,x))dxdydzdw. 
JS 2 xS 2 xS 2 xS 2 

Now recall the reproducing kernel formula 

P l ({x,y))P l {{y,z))dy = 1 ^-P l {{x,y)), 
whence, 



cum 4 (h l . 2 )=48\-^-\ f P?((w,x))dxdw = 48 ^\ 
121 + 1) Js 2 xS 2 (2i + iy 

We then have 

*v(Vfa w ,z) = o( J ° m '^, )=o(r'), 

y{^ar(/i i;2 )} / 

by the result mentioned above. □ 
The next Proposition covers the case q = 3 and follows easily from results 



m 



[13 E]. 
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Proposition 2.2. As I —> oo, we have 



Proof. Note first that 

^;3 = / ff{x)dx = y ai mi aira 2 ai m3 I Yi mi Yi m2 Yi m3 dx 
Js 2 m f^ m „ Js 2 



777177127773 



Illl 



(7713 

77117712 7713 

I I I \ 1 



Z Z I \ f I I I \ /(2/ + 1) 3 



mi m2 7713 / I Om 4-7T 



where 



; V^f' 

fill 



a/(21 + I)" 3 V m l m 2 ™3 



771177427713 

is the normalized bispectrum (see [21 fTH| HH] for definitions and further dis- 
cussions); the symbol on the right-hand side denotes the so-called Wigner 
coefficients, for which we refer to [23] and the Appendix below. It was shown 
in [13] that 

EI m = , Elf u = 6 and ^ d N(0, 6) as I ->• 00 . 



Hence, in view of Lemma IA.ll 



3 l2 / I I I V 3 2 V3 

x 



lim l 2 Ehf. 3 = hm Z 2 f ! ! ! Y 



l^oo l ' 6 2vr l^oo \ J 2vr ttv 7 ^ tt 2 ' 
as claimed. □ 

For the higher order terms, we have the following: 

Lemma 2.3. For all q > 3, as I —> 00 we /lave 

(5) =q\f Pl((x,y)) q dxdy = O . 

The proof of Lemma [2.3l is given in Appendix [B] for completeness. In fact, 
it is possible to obtain the asymptotics of the LHS of ([5]): 



,4, j.. . lo 8 l 



/ Pi((x,y))*dxdy ~ c 4 • . 



and for q = 3 or g > 5 we have 



Pi((x,y)) q dxdy ~ Cf/ . _ 
S 2 xS 2 ' 

for some constants c q (in particular, we get rid of the logarithm on the RHS 
of (jSJ) ) . We use the precise estimates in our forthcoming paper [29] in order 
to compute the variance of {D{\. 
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3. On the proofs of the main results 

Our aim in this Section is to discuss heuristically the expansion of the em- 
pirical measure and process (H])-© into Hermite polynomials. Recall first 
that {fi(x)} is defined pointwise as a Gaussian variable, whence the func- 
tion I(fi(x) < z) belongs for each x,z to the L 2 space of square integrable 
functions of Gaussian variables. In particular, we have that 

I(Mx)<z) = J2 J T-H q (f l (x)), 
where the right hand side converges in the 1? sense (i.e. 

2 



lim E < 



Y,*@-H q Wx))\ =0, 

q=Q q ' 



uniformly w.r.t. z,x). It is possible to provide analytic expressions for the 
coefficients {J q (.)} , indeed for q > 1 

J q {z) = j \u < z)H q {u)4>{u)du = f ^^(u)^ ^(u^^du 

= {-l) q f |^(n)|dn=(-l)^W( z ) = (-l)^-i) (z)) 

where <p{z) := d<fr(z)/dz is the standard Gaussian density. Hence for instance 
JoO) = Ji(z) = -</){z), J 2 (z) = -z(f>(z), J 3 (z) = (1 - z 2 )(/)(z) and in 

general 

(6) J q (z) = -H q (z)cf>(z). 

For z = it follows that Ji(0) = -^vr)" 1 / 2 , 



(-1)" 



■W = 
We have, formally 

I(Mx) < z)dx 

s 2 



x ,Jq-2)l\ q = 2k + l,k=l,2, 
q = 2k, k = 1,2, 



/ jt^H q ( fl (x))dx 
47rxd>(z)+ / jr^±H q (Mx))dx, 



whence 



Gi(z) = Vl [ f2^H q (Mx))dx, 
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and heuristically 

Gi(z) = -^Vl [ H 2 {f l {x))dx 
2 Js 2 

(7) + ( '~y W Vt J g2 HM(x))dx + ... 

m = ^ + !Lz*m ^ lhl , _ c z*m vi hl + .... 

ZD 24 

where are as earlier in ([3]). The meaning of the exchange between the 
integral and the series is discussed below. From the results in the previous 
section, we know that the variance of Vlhi ]2 is given by 

and the variances of the other terms are bounded by 

Var [Vlh^] =o(j), Var [Vlh^] = O f ^ , q > 4. 

In fact, one can compute the precise asymptotic shape of the latter (see the 
discussion after Lemma l2.3[) . 

Our idea in the proofs to follow is hence to show that the term |v / ^iZ;2j' 
is asymptotically dominant in the expansion of Gi, all the other summands 
being uniformly of smaller order. Note that, as remarked in the Introduc- 
tion, the first term in the expansion ([8]), (q = 2) vanishes for z = 0; as a 
consequence, in the Defect statistic D\ the quadratic term 2~ 1 zcf)(z)\/lhi.2 
disappears. From the arguments below it is easy to realize that D\ has hence 
a lower order asymptotic variance than the empirical process; no single term 
is asymptotically dominant and the determination of the limiting behaviour 
becomes much more challenging. The asymptotic behaviour of the sequence 
{D[} is hence delayed to the forthcoming paper [29] . 

4. Proof of Theorem 11.21 

Proof. The key step in the proof is the following uniform weak reduction 
principle: there exist a constant C such that for all < e < 1, as I — > oo, 

Clog 6 I 

which is the statement of Lemma 14.21 We deduce from this that 
sup Gi(z) - J 2 (z)Vlhi ; 2 = o p (l), 

z 

and the result follows immediately from 

Vlhi-2 ^ d N(0, 1) 

(see [13 E]). □ 

It then remains to prove Lemma 14.21 Our approach will follow closely 
the ideas by [5], which were developed in a rather different context (the 



Pr | sup Gi(z) — J2(z)Vlhi-2dx > e\ 



< 
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empirical process for long range dependent sequences on 
the following notation: 



). We will need 



(10) 


Si(z) 


■= Gi{z) 


- h{z)yflhi^ 




Si(zi,z 2 ) 


= Si{z 2 ) 


- Si(zi) , Zi<Z 2 


and 










Jq(Z2,Zt) 




- J q (zi) , Z 1 < Z 2 



Lemma 4.1. There exist a constant C > such that 
(11) 



E [S^z,)} 2 < ^1 [<&(**) - .. 



/ 

Proof. For notational simplicity and without loss of generality we take z\ = 
— oo, i.e. Si(z\,z 2 ) = Si(z 2 ) = Si(z); the argument in the general case is 
identical. We have trivially 

2 



(12) 



E 



s 2 



I(fi(x)<z)-^(z) + ^H 2 (Mx)) 



dx 



E 



S 2 xS 2 



Kfl(x)<z)-<f>(z) + ^-H 2 (f l (x)) 



I(f l (y)<z)-<S>(z) + ^H 2 (f l (y)) 



alxdy 



L 



E 



S 2 xS 2 



Z0(Z 



l{fl{x)<z)-${z) + ^H 2 {f l {x)) 



Kfliy) <z)- Hz) + Z -^lH 2 (My)) 



dxdy 



Now 
(13) 

E 
= E 



Kfi(x)<z)-*(z) 



^3 qi1 ^ Q2 - 



Z0(Z 



H2{fi{x)) 



KMy) <z)- Hz) 



Z(j)(z) 



H 2 {fi(y)) 



J q2 (z) 

92=3 

■^OO Jqi (z) 



°° T 2 (z) 



9=3 



This result is standard, as X^=3 H qi (fi(x)) belongs to the Hilbert 
space of Gaussian subordinated random variables G(fi), with inner product 
(X,Y) := EXY. Plugging flT3|) into $T2§ and using Lemma we finally 
obtain 

2 



E 



s 2 



I(fi(x)<z)-Hz) + ^H 2 (Mx)) 



dx 



r « J 2 (z) a (™J 2 ( Z )} 



x <; sup / P?{(x,y))dxdy 

q>3 JS 2 xS 2 



O 



log/ 

I 2 
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since 



00 J 2 ,U) 



E 

q=3 



Vax(I(//(x) < z)) = $(z)(l - < oo, 



is the variance of a Bernoulli random variable. 
Lemma 4.2. (Chaining argument) For all < e < 1 



□ 



(14) 

and for e > 1 
(15) 



Prjsnp|^)l>4<? fl ° K ' 



Z e 3 



C 



Pr ^ sup |5,(z)| >s} < y <je" 4/3 + 



log 6 / 

f-7/3 



Proof. For (|14p . we follow again closely [5]. We define 

1 



A(z) = + 



1 



2V27T 



x — 1 <j>(x)dx 



-*exp(-— )I(z<-l) 



1 



2V27T 



z 2 1 
^exp(-y) + exp(--) 



I(-l < z < 1) + 



-zexp(-y) + 4exp(--) 



Clearly A(z) is continuous, lim^-oo A(z) = 0, lim^oo A(z) = 1.483943... = 
A(oo). Also, for $(zi,z 2 ) = <&(z 2 ) - $(zi) = Jo (21, 22) we have 

(16) *(zi, z 2 ) + I J 2 (zi, z 2 )| < A(z 2 ) - A(zi) =: A(z 2 , z x ) , for all z 2 > z x . 

by ([6|). Now define refining partitions 

—00 = zo(k) < z\(k) < .... < z 2 k(k) = 00, 

where 



Zi(k) = inf \ z : A(z) > 
implying that 



A(oo)i 



i = 0,...,2 ft -l, 



k{ Zi {k)) - k( Zl ^(k)) < 



L5 



Define ifc(z) by 



Zi k ( z ){k) < z < z ik{z)+1 (k), 



and note that the sequence {z ife ( 2 )(A;)} is increasing, while {-Zj fc ( z )+i(^)} is 
decreasing; then write, for some integer K = K(l) to be determined later 
(17) 

Si(z) = Si(z io(z) (0),z h(z) (l)) + Si(z h{z) (l),z i2(z) (2)) + ... + Si(z iK ( z )(K),z). 



I(z > 1) 
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The idea is to bound uniformly each term of (|17p by means of (jlip and (|16p . 
For the last term of (|17p we have, by a simple algebraic manipulation 



\Si{ z i K {z){K), z)\ < SlU {I{zi K {K) <h{x) <z iK+l (K))}dx + A7r^(z iK (K),z iK+1 (K)) 

{K))\h l;2 \ 

< v^|^(^ K ( 2 )(K),^( 2 ) + i(^))|+2x47rV/A(^(K),z iK+ i(K)) 
+VlA(z iK (K),z iK+1 (K))\h l . 2 \ 

-8vr x 1.5 



< \Tl \Si(z iK{z) (K), z iK{z)+1 (K))\ 



2 K 



Using the latter together with ()17p . we may bound the latter as 
Prjsup \Si(z)\ > ej < Pr (jSi(z io{z) (0), z il{z) (l))\ + ... + \Si(z iK(z) (K), z 1k{z)+1 (K))\ > |) 



/ nM X 1.5 /rl.5 ., . e 

Pr ' ^ 2 k + 2K \ h &\ > 2 



Since YlT=o £ /(^ + 3) 2 — e /2> we ma y bound 

(18) Prjsup 1^(^)1 >s\ < Prjsup |^(z lo(z) (0), z il(z) (l))| > || 



+ Pr| B? |5 I (^ w (l)^ W (2))| > Yg 
+ Pr ^sup\Si(z iK(z) (K),z iKiz)+1 (K))\ > | 1 

+ Pr jVll||^ ;2 |>£-^Vl|. 

Here, using the idea of refining partitions, we effectively reduced the supre- 
mum over a continuous set (R) to the supremum over the possible nearest 
neighbors of z in these partitions. In view of (jlip . the Chebyshev inequality 
yields 

(19) Prjsup \Si(z ik(z) (k),z ik+l{z) (k + l))\ > ^ ^ 2 

2 fc + 1 -l , 

(20) < J] Prj|^(^ + l)^i+i(A: + l))|>7^)2 
Clog / (fc + 3) 



< 



Z e 2 



i=0 

4 2 J2 1 + 1)) - + 1))] < clogl {k + 3)4 



i=0 



We choose 
(21) 

so that 



K 



48vr 



log 2 



-1.5 



48 



-7T 



+ 1, 



e e 12vr 



<±ovi /- /7- L -<J t t J-^' 1 /7 t 

e 2^ 32vr ' 2 2 X 4 
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Hence we obtain 
(22) 



Pr L/l^ \h. 2 \ > - - < Pr \h l]2 dx\ > ~1 < P*{\h li2 dx\ > 8vr} 

* tk riE i ii/2 ^y = w? riyar { /i/2 ? * w? ri - 

To conclude the proof of flU]), we plug the bounds CED and d22D into (fT5|) ; 
we obtain 



A' 



n ) i • i i / Clog^(A; + 3) 4 V 
Pr^sup|S^)l >4 < — A~ > ^ + 



£ 

fc=0 



57T 



< c\ogi{K + 3) 5 | r 1 



5e 2 



(8tt)< 



< ^l^l + ij^omeCC'.C" >0. 

For the proof of (|15p . we repeat exactly the same argument as before, re- 
placing K in (f2Tj) by 



= log 



487T 

^73 



e 12vr n e e 1 / 3 e , 

tttVL > — > - , because e > 1. 

2 2^ 2 4 ~ 4 ' 



It follows easily that 

Pr \Vl 



< Pr 



< Pr 



rl-5 

w 

£ l/3 



/ {/ 2 (x)-l}cte 
Js 2 



e 12vr /- 



32vr 



{/, 2 (x)-l}dx 



> 



5-' 



{f?(x)-l}dx 



> Sire 2 ? 3 < 



-4/3 



(8tt) 



and 



T> ) |C, •! I / C7lQg/(^ + 3) 5 ^g- 4 / 3 ; _x 

Pr <^ sup > e > < — ; — \-C - — I 



I 



5e 2 



ITT 



< 



C" flog 6 / 



I I e 7 /3 



+ £ 



-4/3 



some C, C,C" > 0. 



□ 



Remark 4.3. The proof above is uses the integral 



\x 2 — l| (p(x)dx 



(x 2 - 1) ^(x)dxl(z < -1) 



+ 



(x 2 - l) <j)(x)dx+- I (l - x 2 ) (j)(x)dx \ I(-l < z < 1) 
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2 , 
1 



— 1 



(x 2 — l) <f>(x)dx- 



(1 - x 2 ) 4>(x)dx^ I (x 2 - 1) 4>(x)dx \ I(z 



2V2vr 



- 2 exp(--)I(z<-l) + 



z 2 1 
^exp(-y) + exp(-- 



I(-l < z < 1) 



+- 



1 



2V2vr 

which follows directly from 



z 2 1 
^exp(-y) + 4exp(--; 



I(z > 1), 



(x 2 - l) <f>(x)dxl(z < -1) = [x(f>(x)]i 



-i 



'l - x 2 ) (f)(x)dxl(z < 1) 



and 



(x 2 - 1) 0(x)dxl(z > 1) 



^exp(-y)I(z<-l), 



I(*<1), 
I(z > 1). 



z 2 1 
^exp(-y) + exp(--^ 



z 2 1 
-zexp(-y) + exp(--) 



,4 /so 



|x 2 — l| (fi(x)dx 



4e-l/2 



2V2tt 

5. Proof of Theorem 11.31 
In view of (|l(jp . we can write 

[I* 



0.48394. 



[Lr] 



= W A L(^;r) + W BL (z;r). 
We shall prove that, as L — > oo, 

WAi(^;r) =5- Woo(z;r), supsup Wb£(z; r) = o p (l). 

2 r 

5.1. Step 1: the proof that Wal{z;t) =>■ Woo(z;r), as L — >• oo. 

Proof. To prove convergence of the finite-dimensional distributions for Wal(z; 
it is enough to note that, for n < r2 < ^3 < r4 

7 f ^ [Lr] 
i*W^(*;r) " J2(Z 



and 



J 2 (. 1 )J 2 (. 2 ) [ g / ^ 2 ^ j2(zi)j2(z2)rii 



z=i 



cum {Wal(^i; n), Wal(z 2 ; r 2 ), Wal(z 3 ; r 3 ), Wal(z4\ r 4 ), } 



L 2 



[in] 



l 2 cum,i{hi 2) — > 0, as L — )• 00. 



i=l 



The multivariate extension is trivial. To establish tightness, in view of the 
results from [3] (see also [18j). it is enough to prove that, for all n < r < 
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r 2, z i < z < z 2 , t = k/L, k = l,...,L, some C > 0, there exist a finite 
measure /i : £>([0, 1] x R) — > R + such that 

E { [W AL (z2;r)-W AL (z 1 ;r 1 )} 2 [W AL (z 2 ;r 2 ) - W AL (zi;r)} 2 } 



E 



J 2 {z 2 ) - JgCgl) 



[Lr] 
£ V^ ;2 
«=[i>l] 



-I 2 



Vl 



[Lr 2 ] 

J2 V^/l,; 2 

«=[Lr]+l 



(23) 



£ { [VFal r 2 ) - W AL (zi ; n )] 2 [W AL (z 2 ;r 2 ) -W AL (z; ri )} 2 } 



E 



J 2 (z) - J 2 (zi) 



[Lr 2 ] 

52 vth h2 

l=[Ln] 



2 r 



J 2 {z 2 ) - J 2 (z) 



[Lr 2 ] 

52 vth, )2 

= [Ln] 



(24) < C{fi([ ri ,r 2 ] x [z l; z 2 })} 2 . 

In our case, both these bounds follow easily considering the measure 

M([ri,r 2 ] x [ Zl ,z 2 ]) := (r 2 - n) {A 2 (z 2 ) - Ai(zi)} , ^(R x R) <2.25 < oo. 
Indeed, for (|23p . it is enough to exploit independence over I to show that 



E 



J2(Z 2 ) ~ J 2 {ZI) 



Vl 



[Lr] 

52 ^lh h2 

l=[L n ] 



-l 2 



^2(^2) - h{zi) 

VI 



< 



2 [ir] 



[Lr 2 ] 

^ ^; 2 
= [Lr]+l 

2 [^ 2 ] 



«=[Lri] J |_ v v J l=[Lr]+l 

< (r - n)(r 2 - r) { J 2 (z 2 ) - J 2 (^)} 2 < (r 2 - n) 2 {A 2 (z 2 ) - Ai(zi)} 2 . 



For (12411 . we have 



E 



J 2 (z) - J 2 {zi) 

Vl 



[Lr*} 

52 Vih 

l=[Lri] 



1:2 



J 2 (z 2 ) - J 2 (z) 

Vl 



[Lr 2 ] 

52 f lh w 

l=[Ln] 



{Mz) - J 2 ( Zl )} 2 {J 2 (z) - J 2 ( Z l )} 
L 2 



2 [Lri] 



+6 



iMz)- Mzi)} 2 {J 2 (z) - J 2 (zi)} 
L 2 



52 E[vih, 2 ] 

l=[Lr x } 

2 l Lr 2} 2 _ 

52 E{Vhh h -2} E{Vhh h - 2 y 



h<h=[ Lr i} 
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2 r r- i 4 



< C uniformly w.r.t. /; note also that L 1 < 
(r 2 — ri), for all ri, r 2 in {1/L, 2/L, 1} , whence 

{Mz) - J 2 ( Zl )} 2 {J 2 (z) - J 2 ( Zl )} 2 g ^j^j 4 



L 2 



< C(r 2 -r 1 ) 2 {A 2 (z 2 )-A 1 (z 1 )K 



/=[Lri] 



»]-««»■{*(.)-«■ [ g £ { Ak , li2 } 2 £ {^ i2 } ; 



L 2 

< (r 2 -r 1 ) 2 {A 2 (z 2 )-A 1 (z 1 )} 2 , 
as needed. 



h<h = [L ri ] 



□ 



5.2. Step 2: the proof that sup z sup r Wbl(z; r) = o p (l) , as L -> oo 

Proof. For fixed z, r it is enough to show that 

sup sup |Wbl(;z; r)| = o p (l) , as L — > oo. 

0<r<l — oo<z<oo 



To establish this result, we note that 



El sup sup \W BL (z;r)\} < -^E< 



0<r<l — oo<2<oo 



Vl 



[Lr] 



sup ^ SU P \ S l( Z )\ 



0<r<l — oo<z<oo 



< 4fE^{ su P is,(z)i) 

V I— oo<z<oo J 



C A log 6 ? „log 7 L 



c- 



— >■ , as L — > oo 



because 



£<! sup \Si(z) 

-00<2<00 



< 1 + 



< C" 



Clog/ 



,log^ 
I 



-4/3 , lQ g^ 



/l [ir] \ 

Wbl( z 'j r ) = O p ([EW 2 BL (z-r)] 1/2 )=O p UY, ES i( z n 



• 1 I Lr l i / 



Op(l). 



□ 
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Appendix A. Background on Wigner Coefficients 



Throughout this paper, we made a heavy use of Wigner's 3j coefficients. 
In this appendix, we review briefly some of their features and provide some 
results on their asymptotic properties. We refer to [25], [23] and [I] for a 
much more detailed discussion, in particular concerning the relationships 
with the quantum theory of angular momentum and group representation 
properties of 50(3). 

We start from the analytic expression (valid for mi + m 2 + m 3 = 0, see 
j, expression (8.2.1.5)) 



h h h 
mi m 2 ^3 



(-i)' i+mi V^3+T 



(h+h-kWi-h+hWi-h+ky. 



-11/2 



(h+h+h + iy. 

(l 3 + m 3 )l(l 3 -m 3 y. i 1/2 



(l 1 +mi)\{h-m 1 y.{h + m 2 y.{h 



m 2 ) 



(-l) z (h + h + mi- z)\(h - mi + z)\ 
zl(h +h-h~ z)Kh + m 3 - Z y.(h -l 2 -m 3 + z)\ ' 



where the summation runs over all z's such that the factorials are non- 
negative. This expression becomes much neater for mi = m 2 = m 3 = 0, 
where we have 

^1 h h 




(25) 



1) 



s l+'2-'3 



, for h+l 2 + k odd 
[(h+h+h)/2]\ / (h 

i(ii+i2-h)/m(h-i2+h)/m(-h+i 2 +h)/2\\ \ 

for li + l 2 + l 3 even 



1/2 



Some of the properties to follow become neater when expressed in terms 
of the so-called Clebsch-Gordan coefficients, which are defined by by the 
identities (see [23], Chapter 8) 



(26i 



h 



h 



mi m 2 —m 3 
(27) CgE 



1 



hm2 



(-l)' 1 -' 2+m V2/3 + 1 



Qhmz 



V2/ 3 + 1 h-mih-mi 

h h h 
mi m 2 —m 3 



We have the following orthonormality conditions: 



(28) 
(29) 



{ ^ / l 1 m 1 l2m2 Lj hmil2m2 

mi, ma 



Hi rmf 
°l °m > 



fUlj z m 2 
"mi 0rri2 • 



Now recall the general formula ([24]. eqs. (5.6.2.12-13)) 



5-' 



Y hmi {x)...Y lnmn {x)dx 
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47T 

21 -\- 1 / / / / | l\m\hm2^ L\M\l-3,mz , ' ,K - J L n -zM n -zl 



— ^2 XT \c^ lM } Cr 2 ^ 2 , 



L\...L n -3 Mi...M n -s 



(47T) 



n-1 1 L 'h0h0 Ly L 1 0h0--- L 'L n -3Qln-i0 



Hence we have 



i P W )dt = V WTW I s S° {x)dx = ^ { c ttoC L L f oh o---C[ 

V ' Li...L n -3 

_ 1 /^LiLa-.-iO \ 2 

_ 2Z + 1 ^ \ ho« 2 o...«„-i0/ ' 

Li...L n _3 

in the notation of [15], [16]. Special cases are 

f P?(t)dt = \ ( P? '{cos $)d cos $ 
Jo 2 j 



4tt 



2-7T /"7T 

Y^(tf,^) sin tfdtfc^ 



o Jo 



V ( 2/ + 1) 

1 l>io \ 2 _ p i i V 

" 27+1 I ioio J " V J 



Pl{t)dt — {Qo1o;o} - q| , 1 )^ {QoFqClozo} 

L=0 L=0 
2/ 

= 5^1 E(" + »K» + i) (J 

L=0 V 

(30) = + o ) • 

L=0 V 7 

compare [24J, equation (8.9.4.20). The following result is certainly known, 
but we failed to locate a reference and we provide a proof for completeness. 

Lemma A.l. 1) As I — > oo 

2 



lim / 2 I I I I ) = — = ~ 0.367. 



I I I Y 2 



1-+00 V / 

5j For a// 2 = 1,2, .., and even L = 2, 21 — 2 we have 
I I L\ 2 2 1 



; " 7/i X vr X L(2/-L)V2(2/ + L)V2' 



where 



0.596 = 1.09" 6 < 7, x < 1.09 5 = 1.539. 



3) For L = 0,21 we have 

1 1 °V = _JL_ (' z 2/ V = ^ u + or 1 )! 

J 22 + 1 ' V ; ^(4l + l)Vl X ] ' 
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Proof. 1) Note that (see ([25]) and [24J, equation (8.5.2.32)) 



I I I 




(-iy/ 2 (3Z/2)! 



[II 



(3Z + 1)! 



[(Z/2)!] 3 

Hence, recalling Stirling's formula 

l\ = V2^(l) l+1/2 exp(-Z) + Oir 1 ), 



1/2 



we obtain 

ItJ{1 o o 



lim f 

y/^Y l^oc (1/2) 



(3l/2f l+1 e- 31 

3l+3 P -3l 



l3l+3/2 e -3l 



(31 + l)3/+3/2 e -3Z-l 



\3l+l 



27TI-MO (l/2) 3l + 3 



131+3/2 



- lim 3 3/+1 

7T Z-»oo 



(3Z + l)3«+3/2 
J3Z+3/2 



6e 

, lim 3 

3 I 7T Z— >oo 



3/ 



(3Z + l)3'+3/2 
- Z 3Z 



(3/ + l) 3 ' 



6e 

, lim 

33/27r l— yoo 



(l + l- 1 ) 31 



ttV3 



0.367. 



2) From [23], equation (8.5.2.32) we have 



I I L V _ f (Z + L/2)! 



((L)!) 2 ((2Z-L)!) 

; \((L/2)!) 2 (Z-L/2)!/ (2Z + L + 1)! ' 

We use repeatedly Stirling's formula, for n = 1, 2, ... 

In 1 1 

1 < exp( ) < '-== < exp( ) < 1.09; 

also we write CL n x b n for sequences such that a n /b n ,b n / a n — 0(1) • Hence 
we have 



I I L V _ f (l + L/2)\ 



; -\((L/2)\Y(l-L/2)\ 



(L\) 2 (2l-L)\ 



llL e~ 2l ~ L e - 2L e~ 2l+L 



(2Z + L + 1)! 
(I + L/2) 2i+2i+1 (L) 2i+1 (2Z - L) 2l - L + 1 / 2 



2vr e - 2L e~ 2l + L e -2*-£-i (L/2) 2L + 2 (l - L/2) 2l ~ 2L + l (21 + L + 1)^+3/2 

_ llL xe (1 + L/2) 2l+2L+l (L) 2L+1 (2l - L) 2l - L + 1 ' 2 
~ 2vr (L/2) 2L + 2 (l-L/2) 2l - 2L + l (21 + L + l) 2l + L + 3 / 2 ' 



where 



0.596 = 1.09~ 6 < exp(-— ) < j lL < exp(— ) < 1.09 5 = 1.539. 
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(I + L/2) 



21+2L+1 



(L) 2L + 1 (2l-L) 2l ~ L + 1 / 2 



{L/2) 2L + 2 (l - L/2) 2 '- 2i + 1 (21 + L + l)«+i+3/2 

2 2L+2 2 2Z-L+l/2 y + L ^2l+L+1 tfL+ly _ L ^2l-L+l/2 

2 2i + L + 3 / 2 L 2L + 2 (Z - L/2) 2l - L + l (I + L/2 + i/2)2«+i+3/2 

2 1 26" 1 



L(Z - L/2)V 2 (/ + L/2) 1 / 2 (i + - u J m )2i+L+3/2 L{ i _ L /2)l/2(/ + L/2 )l/2 ' 
whence the proof of the first is completed. 

3) The first part is equation (8.5.1.1) from (|24j). For the second part, it is 
sufficient to recall from [23] equation (8.5.2.33) to deduce that 

2 



I I 21 




1 



(20! 

41 + 1 I IW. 



mm 

(4/)! 



1/2' 



1 



(20 



21+1/2 



41 + 1 



2^+1/2^+1/2 



2Z+l/2/o/\2H-l/2 



2vr(20 2 ' +1/2 (2/) 
(4/)«+ 1 /2 



1/2- 



V2 



as claimed. 

Remark A. 2. Lemma \A . 1\ implies immediately 



□ 



/ I L 




< 



1 



Also, for L = 1 we obtain 



L 2 {4l 2 - L 2 } 
0(r 3 ) , forL = 21. 

1 



, for L = 2,4.... ,21 -2 



L(l-L/2) l / 2 (l + L/2y/ 2 i(i/2) 1 /2(3i/2)V2 ^2' 
leading to the special case 



lim 

Z— »oo 



2 




_vrV3/ 2 _ 





I z Z 




see also |24| . equation (8.9.4-20) for different asymptotic approximations. 



Appendix B. Proof of Lemma 12731 

Proof. First, we note that, since for every t S [—1,1], |-F$(t)| < 1, it is 
sufficient to prove the statement of Lemma 12.31 for q = 4. In this case, we 
have 



rl 

E[hf. 4 ]= 41 Ff(t)dt = 12/ if(cosi?)dcostf 
Jo Jo 



12 

47T 



27T /*7T 



J 



P ; 4 (cos d)d cos 



12 



4tt 



4vr V I 2Z + 1 



Y it{x)dx. 
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Now from (|30j) and Lemma lA.ll 
/ / P; 4 (cos7?)dcosi?# = 47r / i?(t)dt = V(2L + 1) ' * 



L=0 



2L-2 



£ pTHF + g (2i + 1) l»{4P 1 -l'} + 0(rS> - 

For L < / 

E^z^^^E^^, 

L=2 L J L=2 

while for / < L < 21 

+ } L*(2l + L){2l-L} ~ 2P ^ {21 -L} " Z2 • 

Note also that 

2Z-2 2/-2 
L=2 L J L=2 

so this order cannot be improved. Also, using the second part of Lemma 
IA.1I we have easily 

i 



< 4 ' + 1 >(o o o) = ^pW 1 + 0(rl)} = 0(r2) ' 



□ 
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